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of ag(Mz) is in good agreement with the current world average. 
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1 Introduction 



Precise data on e + e~ annihilation into hadrons have provided detailed experimental tests 
of QCD and one of the best opportunities to measure the strong coupling constant as- A 
particularly well suited observable is the energy-energy correlation (EEC) pQ, defined as 
an energy-weighted correlation for the cross section corresponding to the process e + e~ — > 
h a + h b + X as 

— , ^ = — l -7^T da e+e-~>h a h b +x5(cOSX + COs6 ab ) , (1) 

a T d cos x °t *-? J Q 2 

a,b 

where E a and E b are the energies of the particles, Q is the centre-of-mass energy, 9 a i> = n — x 
is the angle between the two hadrons, and ot is the total cross section for e + e~ — > hadrons. 

The two-hadron cross section e + e~ — > h a + h b + X depends on the fragmentation 
functions of the partons into the final-state hadrons. However, thanks to the momentum 
sum rule 

1 

dxxD h/q (x,fi 2 F ) = 1 , (2) 



o 



EEC becomes independent of them, and can thus be computed in QCD perturbation theory. 

Theoretical calculations 01 H] for the EEC function have been performed up to next- 
to- leading order (NLO) accuracy in QCD [5] allowing a detailed comparison with the 
available data. 

As is well known, fixed-order calculations have a limited kinematical range of applica- 
bility. In the back-to-back region, defined by 9 a b — ► n (x —> 0), the multiple emission of soft 
and collinear gluons gives rise to large logarithmic contributions of the form a<? log 2 ™ -1 y, 
where y = sin 2 x/2. As y decreases, the logarithms become large and therefore invalidate 
the use of the fixed-order perturbative expansion. 

These logarithmic contributions can be resummed to all orders • The resumma- 

tion formalism is very close to the one developed for the transverse-momentum distribution 
of high-mass systems in hadronic collisions When the transverse momentum of the de- 
tected final state is much smaller than its invariant mass Q 2 , large logarithmic contributions 
ag log 2n_1 qr/Q 2 arise which must be resummed to all orders. 

The coefficients that control the resummation at a given order can be computed if 
an analytical calculation at the same order exists. In the case of hadronic collisions, the 
complete form of the logarithmically enhanced contributions has been computed |15j . In 
this paper we present the result of a similar calculation, performed for EEC. Our calculation 
allows us to fix the still missing coefficients at C?(a|) and to extend the resummation for 
this observable to full next-to-next-to-leading logarithmic (NNLL) accuracy. We also study 
the numerical impact of our results and present a comparison with LEP and SLC data. 

The paper is organized as follows. In Sect. 2 we review the resummation formalism and 
we discuss the results of our calculation. In Sect. 3 we present numerical results, and we 

TThe role of the transverse momentum is played, in the case of EEC, by the variable q\ = Q 2 sin 2 (x/2). 



1 



also consider the inclusion of hadronization effects. In particular, we perform a fit to OPAL 
and SLD data. 



2 Soft-gluon resummation 



The EEC function can be decomposed as 



1 cffi 1 d£( ros -) 1 d£( fin -) 

+ — ( 3 ) 



cry d cos x &T d cos x d cos x 
the first term on the right-hand side of Eq. Q contains all the logarithmically enhanced 
contributions, ag/?/log m ?/ at small y and has to be evaluated by resumming them to all 
orders. The second is free of such contributions and can be computed by using fixed-order 
perturbation theory. 

The resummed component can be written as [T31 ITl] 

i ^vCrcs.) f)2 roc 

— -, = ^H(a s (Q 2 )) / dbbJ (bq T )S(Q,b). (4) 

ot d cos x o Jo 

The large logarithmic corrections are exponentiated in the Sudakov form factor 

2 d£ L,_ ^ ,_Q 2 



S(Q, b) = exp { - 



A(a s (q 2 )) ki^ + B(a s (q 2 )) 



(5) 



The Bessel function Jq(&?t) and &o — 2e 7E have a kinematical origin. 



The resummation formula in Eq. (J1J) has a simple physical interpretation. When the 
triggered partons are back to back, the emission of accompanying radiation is strongly 
inhibited and only soft and collinear partons can be radiated. The function H(as(Q 2 )) 
embodies hard contributions from virtual corrections at scale q ~ Q. The form factor 
S(Q, b) contains virtual and real contributions from soft (the function A) and flavour- 
inclusive collinear (the function B) radiation at scales 1/b^q^Q. At extremely low scales, 
q&l/b, real and virtual corrections cancel because EEC is infrared safe. 

The functions A, B and H in Eqs. (|4|5|) are free of logarithmic corrections and can be 
computed using a perturbative expansions in as: 

oo 

as 

71 



A ^ = E (?) aH ' (6) 

n=l 

oo 

^ = £(v) v,) • (7) 

n=l 

oo 

H(a s ) = i + £(«)V"> . (8) 



71=1 



By explicitly performing the q 2 integration in Eq. (|SJ) the form factor can be recast in 
the following form [THl H3 [TH] : 

S(Q,b) = exp{L g^asPoL) + g 2 (a s poL) + a s g 3 (as(3oL)...} , (9) 
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where as = as/ 71 " an d the large logarithm L = \ogQ 2 b 2 jb\ at large b corresponds to the 
logy, which becomes large at small y (the limit j/ < 1 (qx <C Q) corresponds to Qb ^> 1 
through a Fourier transform). 



<?i(A) 



The explicit expressions of the gt functions are*: 
A« A + log(l - A) 



A 



Po 

2 (A) = -5-log(l-A 



+ 



<? 3 (A) 



A) 
Pi 

A(3) 



A 



1- A 

log(l - A) 
1 - A 



+ log(l-A) + 



AW 
Po 



1-A 



- I \ )log% 



A 



1-A 



A 2 



+ 



2/?o 2 (1 " A) 2 
/ A 



\ A^fo fA(3A-2) (1 - 2A) log(l — A) 



A 2 



1-A 



+ 



A) 1-A 

log(l ~ A) 
1-A 



ft V2(l-A) 2 



(1-A) 2 



A 2 



2 (1-A) 2 



, 2Q 2 
/4 



A 2 



%l V 1"A + A (1-A) 2 



A 



+ ^ (1) ^r 



+ 



2A 



(1-A) 



log(l - A) 



A 2 1-2A 
2Af(T^Ap 



log 2 (l-A) + log(l-A) 



. Po 



+ 



A 4 (i - A) 



A 



:(/3 A(2-3A)+A 2 A) 



2A 4 (1-A) 2 ' 

and the coefficients of the QCD (3 function are defined as: 



Po = - (11C A - 2n/) , 

1 Z2857 3 1415 
^ 2 " 64 " ~^4~ 



A = i (17C 2 - hC A n f - 3C F n f ) 



C\n f - ^C A C F n f 4 C 2 F n f 4 ^C A n 2 f + yC>7iJJ 



(10) 



(11) 



The functions g%, g2, g% control the LL, NLL, NNLL contributions, respectively. The 
coefficients and B^ 1 ' were computed a long time ago ^3] and are the same as 

appear in the quark form factor in the transverse momentum distributions in hadronic 
collisions. They read 

A® = C F 



fl« = -°- Cf 
2 



/67 



A^ = -(C A (--^ 
2V A \18 6 



(12) 



By using for a F the NLO expression 

47ra 2 



NLO 

(Jrp 



Q- 



(13) 



■^Throughout the paper we use the MS renormalization scheme. 
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the coefficient is [1 



11 



TP 



c F — + — 



(14) 



The coefficient has been obtained recently, as the leading soft term in the three-loop 
splitting functions [T§1 |2U]: 



A (s) 



C 2 



245 67 



11 



11 



55 
24 



+ 2C; 



+C A n f 



209 10, 7, 

1 C2 C3 

108 9 s 3 s 



1 2 
27 f 



(15) 



where ( n is the Riemann ( function (£ 2 = 7r 2 /6, (3 = 1-202. 



In order to be able to perform the resummation up to NNLL accuracy, only the coeffi- 
cient B^ is lacking. There have been in the past several attempts to obtain a numerical 
value for this coefficient ^HIIS]- I n this work we will present the analytical result for B^ . 

A direct way of extracting the resummation coefficients consists in comparing the loga- 
rithmic structure of a fixed-order perturbative calculation of the EEC, with the expansion 
of the resummed formula in Eq. (jlj). The expansion up to 0(a\) reads 



1 d£( res -) 
ot d cos x 



J_ fas 



-A^\ogy + B^] + {^-) 
+(3 AW) log 2 y 



+ (-A^ - /3 BU + (5(D) 2 - A^H^A log;, 



)\ 2 



+ 0(4) 



(16) 



where we have set fiR = Q. 



An analytic calculation of EEC at NLO (i.e. up to 0(a\)) would allow the extraction 
of the coefficients A^\ B^\ H^ 1 ', A^ and B^. However the full analytic result is not 
really necessary to this purpose: it is sufficient to compute its small-?/ behaviour. 

The strategy to obtain the small-?/ behaviour is the one applied for a similar calculation 
in the case of the transverse- momentum distribution in hadronic collisions ^5] • The singular 
behaviour at small y (qx) is dictated by the infrared (soft and collinear) structure of the 
relevant QCD matrix elements. At Q(ais) this structure has been known for a long time 
[22] • In recent years, the universal functions that control the soft and collinear singularities 
of tree-level and one-loop QCD amplitudes at 0(a>§) have been computed [23 121] • By using 
this knowledge, and exploiting the simple kinematics of the leading-order subprocess, we 
were able to construct improved factorization formulae that allow the control of all infrared 
singular regions, avoiding problems of double counting [T3]. We have used these improved 
formulae to approximate the relevant matrix elements and compute the small-?/ behaviour 
of EEC in a simpler manner. 
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Compared with the calculation of Ref. ^3], i n the case of EEC there is an additional 
complication. The definition of EEC in Eq. (0) implies that a sum over all possible par- 
ton pairs has to be performed. Thus an infrared-finite result can be recovered only after 
summing over all the correlations. 

For hadron-initiated processes, the coefficient Ba is generally dependent on the re- 
summation scheme, and on the way the resummation formula is actually organized |18j . 
However, despite these ambiguities, it always has the form [T3] 

B^ F = -\^ + \(3oA F a , a = g ,0, (17) 

where 7a is the coefficient of the 5(1 — z) term in the two- loop splitting function [23 
I2H]. The second term in Eq. (JT7J) depends on the virtual correction to the process aa — > 
F(qx, Q 2 )- Considering the similarity between EEC and the transverse momentum spectra 
in hadronic collisions, it is natural to expect a similar form for the coefficient B^ for the 
EEC, modulo possible crossing effects. 

More precisely, since the leading-order subprocess which is relevant here is the produc- 
tion of a qq pair, we expect 

5 (2) = —1? + C F f3 X. (18) 

Assuming Eq. (JIBJ), a calculation of one of the two colour factors CpT R or CpCa is sufficient 
to fix the coefficient X in Eq. (fTHj) . We have computed both the CpT R and the CfCa 
contributions to Eq. ()16|) and found complete agreement with all known results. Our 
results are also consistent with Eq. (|TH|) and allow us to fix 

5 (2) = -^ 2) +^o(^ 2 - 2 ) > ( 19 ) 

the coefficient 'y^ being 

7f » = Cl (I - £ + 6C3 ) + C F C A g + i£ - 3 C3 ) - C F n,T R g + . (20) 

In principle, since the contribution to Eq. ()16|) proportional to the colour factor C\ has 
not been computed, the result in Eq. (fTUj) is not fully established. However, besides the 
parallel with transverse momentum distributions in hadronic collisions, which strongly 
suggests Eq. (|18j) . there are two additional arguments that confirm it § . The first one 
relies on the correspondence that should exist between our coefficient B^ and the quark 
coefficient in the non-singlet (NS) scheme jTHj: 

*8* = -?? + ™$-\), (21) 

which is expected to directly measure the intensity of collinear radiation from quarks at 
0(0:5). We find that the difference between Eqs. ()19|) and (|2*T|) can indeed be explained as 

§We also note that our result agrees with the one guessed by K. Clay and S.D. Ellis [H] based on the 
similarity of EEC to Drell-Yan and, as far as the Cftif part is concerned, with an independent calculation 

EH 
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a pure crossing effect due to an additional factor present in the phase space in the case of 
EEC. 

Finally, the numerical value of for rip = 5, B^ 2 > = 11.2, is in good agreement with 
the estimate of Ref. [21], B^ ~ 10.7 ' obtained with the numerical program EVENT2 |llj . 

The resummed component obtained in Eq. has to be properly matched to the fixed- 
order result valid at large y. The matching is performed as follows: 

1 c/£( fin -) 
<7t d cos x 

The first term on the right-hand side of Eq. ()22j) is the usual perturbative contribution; it 
is computed with the numerical program of Ref. ^T] at a given fixed order (LO or NLO) in 
ag. The second term is obtained by using the expansion of the resummed component (see 
Eq. ()16|)) to the same fixed order in as- This procedure guarantees that the right-hand side 
of Eq. © contains the full information on the perturbative calculation plus resummation 
of the logarithmically enhanced contributions to all orders. 

As we will show in the next section, our result in Eq. (fTU|) allows us to perform an 
excellent matching between the resummed and perturbative NLO result. 

We finally note that the functions gi are singular as A — > 1. The singular behaviour is 
related to the presence of the Landau pole in the QCD running coupling. To properly define 
the b integration, a prescription to deal with these singularities has to be introduced. Here, 
analogously to what was done in Ref. j2H], we follow Ref. [2H| and deform the integration 
contour to the complex 6-space. 



1 rfS 




1 d^ Tes -> 


ot d cos x 


f.o. 


a T d cos x 



(22) 



3 Phenomenological results 

In the following we present quantitative results at NLL+LO and NNLL+NLO accuracy. 
At NLL+LO the resummed component in Eq. (0J) is evaluated by including the functions 
gi and #2 in Eq. (fiT)|) and the coefficient in Eq. (fTlj) . The finite component in Eq. (}2"2*|) 
is instead evaluated at LO and the one-loop expression for as is used. At NNLL+NLO we 
include also the function g 3 in the resummed component and we evaluate the finite part at 
NLO, with as at two-loop level. 

The NLL+LO results are shown in Figs. ^ and 01 They are obtained by fixing Aqq^ 5 = 
0.1665 GeV, corresponding to as(M^) = 0.130. In Fig. [T]we show the results for = Q. 
The dotted line is the LO result, which diverges to +oo as x 0- The solid line is the 
matched result, and the dashed line gives the matching term in Eq. (}2"2*|) . Note that we 
plot l/cr T dTj/dx, so that the matching term is actually Y(x) = l/°r dE^ 6n ^/dx- As can 
be observed, the matching term is well behaved up to very small values of x an d becomes 
dominant at larger x, where the fixed-order contribution is expected to control the matched 
calculation. 

^More precisely, the numerical estimate of Ref. [2] is for the coefficient G21, which is related to 
by = \G 21 + JLC F n F - \ (f - C F C A - C(3)C|. 
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The scale dependence at NLL+LO is studied in Fig. where the results for the scales 
H>r = Q/2, Q, 2Q are shown. The lower plot shows the detail of the region 30° < x < 60°. 

The NNLL+NLO results are shown in Figs. 01 and 0] They are obtained using Aqq^ 5 = 
0.23 GeV, corresponding to a$(Mz) = 0.118. As before, Fig. |3] shows the results for 
Hn = Q. The dotted line is the NLO contribution, which diverges to — oo as \ — > 0. The 
solid line is the matched result, and the dashed line gives the matching term in Eq. (|22jl. 
We see that this matching term displays a very smooth behaviour as \ — > 0, and this is 
a further confirmation of the validity of our result in Eq. (|T^|) . The NNLL+NLO result is 
plotted in Fig. H|for = Q/2, Q, 2Q. As in Fig. |21 the lower panel shows the detail of the 
region 30° < \ < 60°. 

We see that scale variations act differently in the low-% and medium-x regions. In 
the region of the peak, lowering (increasing) fipt has the effect of increasing (lowering) 
Q! s(A i i?) an d thus increasing (damping) the Sudakov suppression. This results in the fact 
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that, at NNLL+NLO, the curve at fin = 2Q is higher than the one at = Q/2. This 
behaviour changes as x increases and at x — 60° the curves are in the usual order. At 
NLL+LO these two distinct features (Sudakov suppression and perturbative increase with 
as) are less evident, and thus the scale dependence appears smaller than at NNLL+NLO. 
At NNLL+NLO the scale dependence is about ±8% at the peak and about ±5% at x — 60°, 
giving an idea of the theoretical uncertainty in the resummed calculation. 

We find that the NNLL effect is dominated by the contribution of B^ in the function 
#3. By keeping only the term proportional to in the function g% in Eq. (jlUj) . the 
difference with respect to the full NNLL+NLO result is smaller than 1%. 

Figure 03 shows the NNLL+NLO matched result in the full range of x- We see that, 
contrary to what happens in other approaches to 6-space resummation (221; there are no 
oscillations in the medium-high x region, where the matched result follows the NLO fixed 
order calculation. 
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Figure 5: NNLL+NLO results in the full \ range 



Before moving to the comparison with the experimental data, we want to study the 
convergence of the resummed expression, as a check of the validity of the prescription 
introduced in [2HJ and used to deal with the Landau pole. In Fig. |H1 we compare the 
purely resummed result in Eq. (JU) at NNLL accuracy (that is with the functions g±, g 2 , 
#3 and the coefficient H^ l > included) with its expansion up to 0(a$). As can be observed, 
the expansion converges very rapidly to the resummed result in the region of medium x, 
confirming the validity of the prescription where the fixed-order result dominates (see the 
lower plot for a detailed comparison). Nonetheless, even though the higher the expansion 
the better the agreement with the resummed result at smaller values of \i f° r X ^ 10° the 
fixed-order result is bound to fail, no matter how many orders in perturbation theory are 
included, thus requiring the resummation to all orders. 

We can now perform a comparison of the most accurate theoretical NNLL+NLO results 
with the precise OPAL [3U] and SLD [3T] data. As the extraction of the strong coupling 
constant as is one of the main motivations for the measurement of event-shape observables, 
we perform a fit of the experimental data on EEC leaving Aqcd as a free parameter. We do 
not attempt to produce the most accurate extraction of a$(Mz), since we cannot properly 
take into account correlations between the data points and therefore just add systematic 
and statistical errors in quadrature. 

For the moment we neglect hadronization effects, but we will come back to this point 
below. The reader should keep in mind that the results obtained without including those 
effects should be considered with care. 

In a first sample we include data in the range 15° < x < 120° and fix the renormalization 
scale to fifj ; = Q = Mz- The upper limit is chosen so as to cut the large angle region where 
another resummation would be required. 

The quality of the fit is poor, as can be determined by the value of x 2 /d.o.f. = 5.17, 
with a rather large value for the coupling constant as(Mz) = 0.133 ± 0.002, in agreement 
with similar results found by OPAL [3U| . The uncertainty is dominated by missing higher- 
order contributions, estimated by repeating the fit with ji R = 1/2 {2)M Z - Better fit results 
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are actually found when the renormalization scale is also varied, with a reduction of a 
factor of 2 in \ 2 when /j,r ~ Mz/2 and a slightly lower preferred value of the coupling 
constant. In a second attempt, we include data in the range 0° < x < 63°, to isolate 
the region where the effects of the resummation are considerably more significant. Even 
with the scale fixed to /i« = M z , a very reasonable value of x 2 /d.o.f. = 1.67 is found, 
corresponding to as(Mz) = 0.131 ±0.002. The error is again dominated by scale variations 
fiR = l/2Mz,2Mz- The nice comparison between the resummed calculation and the data 
in the low-angle region is shown in Fig. [7| It is clear that our NNLL+NLO result for the 
EEC can reproduce very well the data up to the lowest measured angle. 

Since the low-x region is particularly sensitive to non-perturbative (NP) effects, whereas 
in the large angle region we may expect non negligible higher-order (NNLO) contributions, 
we have repeated the fit in the range 15° < x < 63°. The value of x 2 /d.o.f. goes down to 
0.66 but the result for cts(Mz) does not change significantly. 

We have also investigated the possible effect of the unknown second-order coefficient 
(see Eq. (jHJ)), by letting it vary together with ces(M z ) in the range 0° < % < 63°. The 
results show that data still prefer a high as(Mz) and a relatively small H^ 2 \ 

Up to now we have considered only the perturbative contribution in the theoretical 
calculation. However, NP contributions are expected to be relevant, particularly for small 
angles [TJ HH1 EH] ■ Thus, following Ref. [21], we include NP effects by supplementing the 
Sudakov form factor in Eq. (JH]) with a correction of the form 

S NP = e -^ aib \l - 2a 2 b) . (23) 
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Figure 7: NNLL+NLO: purely PT fit in the low-x region 

We have performed a three-parameter preliminary fit to the data still in the range 0° < \ < 
63°. We find that the data prefer very small values of the NP coefficient ci2, [ eta | ^0.002. 
We have thus set a<i = to perform the fit. We obtain for as(Mz) and ai the following 
result: 

a s (M z ) = 0.130±g : gg at = 1.5±g;l GeV 2 (24) 

with x 2 /d.o.f. = 0.99. The error is dominated by scale uncertainty, which is estimated as 
above by repeating the fit with /ir = 1/2(2)M^. 

We see that the quality of the fit improves, but the value of as(Mz) still remains high 
with respect to the world average as(Mz) = 0.1182 ± 0.0027 [33J. Moreover, the three- 
parameter fit suggests that in our approach the NP coefficient a\ is more important than 
ci2. We conclude that the parametrization in Eq. (J2"3~|) is not able to fully take into account 
the hadronization effects, particularly at medium and large values of \. The extracted 
"effective" coupling as thus absorbs part of the hadronization effects. We stress that this 
result is not an artefact of the resummation procedure, since a similar effect is observed 
when the large-angle data are compared with the fixed order NLO result. 

A different method to include NP (hadronization) corrections, extensively applied in 
the QCD analysis to event-shape variables at LEP, is to use a parton shower Monte Carlo. 
This method is certainly model-dependent, since different approaches have been developed 
to describe the hadronization process, but has the advantage that the free parameters are 
tuned through fits to large sets of different data. 

In Refs. |3UJ EH] the data for EEC have been corrected at parton level before performing 
the QCD analysis. We have used the parton- level data of Ref. [30J to repeat our fit. The 
hadron-parton correction factors are large, from about 1.5 in the very small \ region to 
~ 0.9 at large \. The quality of the fit in terms of x 2 /d.o.f. is generally worse than 
before, but errors related to the hadronization correction have not been included. The 
uncertainty from hadronization is usually estimated by trying different alternatives for the 
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hadronization correction, and using the spread in the results as an estimate of the ensuing 
error. In Ref. [3U] the hadronization uncertainty on a^{Mz) is estimated to be about ±4%. 

In the kinematical range 15° < \ < 63° we find % 2 /d.o.f. = 3.78 and «s = 0.119±0.001, 
with a very similar result (x 2 /d.o.i. = 5.02 and a s = 0.120 ± 0.001) when the fit range is 
extended to 15° < x < 120°. Finally, when the fit is repeated by allowing the variation 
of the renormalization scale, excellent results are obtained, the best fit corresponding to 
as = 0.1175, x 2 /d.o.f. = 1.36 and /jlr = 0.28 Mz (see Fig. |HJ). Thus, when hadronization 
effects are taken into account using a Monte Carlo, as done in Refs. j301 13T] , the results we 
obtain for a$(Mz) are in complete agreement with the world average. 

We note that the QCD analysis performed by OPAL on the same parton-level data 
gave, instead, ces(Mz) = 0.1321q;qq7 [3Sj, considerably higher than our result. The analysis 

1.0 

0.8 



0.6 

X 

\ 
W 

73 0.4 
b 

\ 

0.2 
0.0 

20 40 60 80 100 120 

0.4 
0.3 
0.2 
0.1 
0.0 

30 35 40 45 50 55 60 

X 

Figure 8: Comparison of NNLL+NLO fit to parton level OPAL data with a$(Mz) = 0.1175 

of Ref. [3U] (as well as the one of Ref. [HI]) used the NLL resummed calculation of Ref. [To] . 
In order to understand the origin of the discrepancy, we have compared our results with 
those of Ref. ^H]- We find that the differences are not negligible, especially in the region 
X^40°, and may explain the above discrepancy. The approach of Ref. [15] was based on 
an approximated analytic evaluation of the 6-space integral of Eq. Q, and suffers from an 
unphysical singularity at very small \- The effect of this singularity may propagate also 
within the range of the fit, thus spoiling the resummed prediction. The other source of 
difference is that the coefficient B^ 2 \ or equivalently G 2 i, was evaluated numerically, thus 
leading to a larger uncertainty in the matching procedure. 





12 



4 Conclusions 



In this paper we have considered the observable known as energy-energy correlation in 
e + e~ collisions. We have provided the complete structure of the 0(a\) logarithmically 
enhanced contributions up to NNLL accuracy, by giving the expression of the unknown 
second-order coefficient B^ 2 \ needed to reach full NNLL precision. 

We have presented perturbative predictions both at NLL+LO and NNLL+NLO accu- 
racy, by showing that the knowledge of the calculated coefficient allows us to perform 
an excellent matching of the resummed and fixed-order calculations. We have studied the 
impact of the results in a detailed comparison to precise LEP and SLC data. A good 
description of the data is obtained but the extracted value of as(Mz) turns out to be high 
when hadronization effects are neglected or parametrized using the form in Eq. (|23|). By 
contrast, using OPAL data corrected at parton level, that were obtained by estimating 
hadronization corrections using a Monte Carlo parton shower, the values of cks(Mz) we 
find are in good agreement with the current world average. 
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